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We consider a system of trimers and monomers on a triangular lattice that describes 
the adsorption problem on the (111) lxl surface of a diamond crystal structure. We 
introduce a mapping to a 3-state vertex model on a square lattice. We treat this adsorption 
problem using the transfer-matrix method combined with a density-matrix algorithm, to 
obtain thermodynamic quantities. 

§1. Introduction 

The two-dimensional adatom eas^ (adsorbate) plays important roles in a wide 
variety of surface phenomena.^' B^tf'tf'til' For example, the growth mode or kinetic 
process of surface dynamics depends on the thermal adatom concentration. Owing 
to recent developments of the surface observation techniques at various length scales, 
quantitative study of the equilibrium adatom concentration and adatom formation 
energy has become possible. The thermal adatom density cq is usually given by 

c = n exp(-E af /k B T), (1-1) 

where no is the density of atomic adsorption sites on the surface, E a f is the adatom 
formation energy, k-g is the Boltzmann constant, and T is the temperature of the 
surface. Using this formula, Tromp and MankosEF have determined the adatom 
formation energy for the Si(OOl) surface to be 0.35 ± 0.05 [eV]. (The "adatom" is 
actually a dimer, as they have also shown.) 

For the Si(lll) high-temperature lxl surface, however, the situation is some- 
what more^omplicated, because the Si adatoms form short-range order of y/3 x 
symmetry. B>' LP 'B>'0> From the half width of the \/3 x y/3 spot obtained by a diffrac- 
tion measurement, the mean size of the ordered island is thought to be several 
nanometers. I12P Since the formula ( |1-1| ) holds only in the low-density regime, where 
the adatoms are almost uncorrelated, it can not be applied to a system with well- 
developed short-range order. Further, on the Si(lll) lxl surface, there are two 
kinds of adsorption sites, called T4 and H3, which are common to the (11 1) 1 x 1 
surface of the diamond crystal structure (Fig. |]). Kohmoto and IchimiyaEP showed 
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Fig. 1. Top view and side view of the (111) surface of diamond structure. 




Fig. 2. Top view of the (111) surface of a diamond crystal structure with adsorbates and the 
corresponding triangular trimer configuration. 

that adatoms are adsorbed at both T4 and H3 sites with a mixing ratio of about 
4:1 j from the analysis of the reflective high energy electron diffraction (RHEED) 
intensity on Si(lll) lxl surface at 900°C. The total coverage of Si is around 0.25 
monolayer (ML) (1 ML = one adatom per lxl lattice point). 

In the present paper, we consider the lxl phase of the (111) surface of diamond 
structure and make a statistical-mechanical study, taking both T4 and H3 sites into 
consideration. In applying our analysis to the Si(lll) surface, we focus on the high- 
temperature lxl phase, and our aim is to determine the relationship between the 
adatom concentration and the adatom formation energylnthis sense, our present 



transition is the main issue. 

As shown below, our problem is basically a trimer-monomer mixture problem on 
a triangular lattice. In the zero monomer density limit, it becomes equivalent to the 
trimer tiling (i.e. close packing) problem, whose exact solution was given recently. 
From the viewpoint of the tiling problem, the inclusion of monomers, i.e., removing 
the close-packing condition, may affect the problem in two ways: (1) the breakdown 
of exact solubility, and (2) the alteration of the phase-transition-related character. 
With regard to (1), we should recall that, in the well-known dimer tiling problem, 
the close-packing condition plays a crucial role for its Pffafian solution. The close- 
packing condition is also essential in the Bethe ansatz solution of_the trimer tiling 
problem. li2P With regard to (2), we should recall a rigorous proof Ej) demonstrating 
the absence of the phase transition for non-zero monomer fugacity in the dimer 
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model. Hence, inclusion of monomers in the tiling problem is a non-trivial and 
interesting problem which deserves detailed study. In the present paper, to obtain 
reliable results for general cases with non-zero monomer density, we use the product 
wave-function renormalization group (PWFRG) method,H which is a variant of the 
density matrix renormalization group (DMRG) method. E-J 1 
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Fig. 3. Three types of elementary vertices. 



2. Generalized trimer tiling problem and its vertex- model 

representation 



We consider a triangular lattice A representing the ideal (non-reconstructed) 
(111) surface of a diamond structure, and we place a dangling bond at each lattice 
point. The adsorption sites are on its dual lattice, A*. We divide A* into two sublat- 
tices, each of which is a triangular lattice. One sublattice represents the adsorption 
sites of T4, and the other represents H3 sites (Fig. [E]). We assume that an adsorbate 
atom on A* must consume three dangling bonds around the adsorption site, and 
that a dangling bond is never linked to two or more adsorbates (the "non-share" 
condition) (Fig. 2). 

We can regard the adatom accompanying three dangling bonds as a "trimer" 
occupying an elementary triangle of the triangular lattice. It should be noted that the 
overlap between trimers is excluded in accordance with the non-share condition. A 
free dangling bond which is not linked to any adsorbate is regarded as a "monomer" 
on the triangular lattice. 

The Hamiltonian of our system is 

n = H( £T 4 ~ 3e s)cT 4 ,i + X!( eH 3 ~ 3e s) c H 3 ,i> (2' 1 ) 

i i 

where the ct 4 ,i (resp., ch 3 ,i) are the lattice gas variables at the i th site of T4 (resp., 
H3). The corresponding adsorption energy is ex 4 (resp., en 3 )- The quantity e s is the 
surface energy per lxl unit cell without adsorbate. The lattice-gas variable takes 
either the value (empty) or 1 (occupied). 

As is shown in Fig. || we introduce three elementary vertices, which we denote 
by A, B, and C. When a T4 site is occupied by an adsorbate (resp., empty), we 
assign (1,1,1) (resp., (0,0,0)) to the legs of the vertex A; the same assignment applies 



4 



N. Akutsu and Y. Akutsu 



(0 0)=1, (0 1)=2, (1 0)=3 

1 2 1 





/ 

o 


\ 







1 


>) 




f 









f 


1 












■X- 




1 1 
1 1 3 



1 3 1 




■X- 



1 3 1 



1 3 1 1 2 

■X" m x 

1 1 1 1 2 

1 1 1 1 

■X- M-X 



1 1 



>P(3£ S - e H ) 



others = 



Fig. 4. Composite vertices and the three corresponding Boltzmann weights, where we make the 
index correspondence (0,0) — > 1, (0, 1) — > 2, and (1,0) — > 3. 



to the vertex B for the H3 site. For the vertex C, there are seven configurations 
corresponding to the state of a dangling bond, as shown in Fig. ||. The seventh 
pattern in Fig. || is regarded as the monomer in the trimer-monomer problem. 




Fig. 5. Temperature dependence of Ct 4 and Ch 3 for a closed packed trimer problem calculated 
using the PWFRG method for eH 3 /eT 4 = —0.5. The exact transition temperature is given by 
—tT i /kB.T m = 4/3 ln(3\/3/4) w 0.34883. The dashed line corresponds to the metastable state. 

From the elementary vertices, we build a composite 4-valent vertex. The associ- 
ated Boltzmann weights define a 3-state vertex model on the square lattice (Fig. 4). 
We thus have a novel mapping of the trimer-monomer mixture problem to a vertex 
model. 

By a standard procedure, 113 we construct the row-to-row transfer matrix T 
for the 3-state vertex model. To obtain_the per-site free energy from the largest 
eigenvalue of T, we employ the PWFRGli^ method to obtain an approximate (but 
fairly accurate) diagonalization of T. Since T is asymmetric, we apply a version of 
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the modified PWFRG method used in Ref. 19), which can deal with an asymmetric 
transfer matrix. 
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Fig. 6. Temperature dependence of concentrations calculated using the PWFRG method in the 
disordered phase. (a) The quantity Ctotai- The ratio r = (3e a — £H 3 )/(3e a — £t 4 ) is 
0.75, 0.5, 0.25, 0, —0.5, — 1, —2, —5, from top to bottom. The bold curve corresponds to the 
case r = 0. (b) Ct 4 (solid curves) and Ch 3 (dashed curves). For Ch 3 (resp., Ct 4 ), 
r = 0.75, 0.5, 0.25, 0, —0.5, — 1, —2, —5, from top (resp., bottom) to bottom (resp., top). The 
bold lines correspond to the case of r = 0. 



§3. Calculated results 

In the limit e s — ► oo, the statistical weight of the monomer decreases to zero. 
Then the system can be described by a close-packed trimer model whose exact solu- 
tion is known.© To see the reliability of our method, we calculate the entropy and 
the sublattice trimer concentrations (in the disordered phase), which we compare 
with the exact solution. H At £t 4 = £h 3 , the PWFRG calculation gives the entropy 
in the disordered phase as 0.0871/cb (per lxl unit cell), which is in good agreement 
with the exact valueB ((l/3)5 sym /fc B = (1/3) ln(3\/3/4) 0.08720). 

In Fig. 5, we give the calculated result of the sublattice trimer concentrations 
(average number of trimers per lxl unit cell) Ct 4 and Ch 3 for eH 3 /eT 4 = —0.5. 
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Fig. 7. Critical curves of the trimer-monomer mixtures calculated using the PWFRG method (solid 
curves). Calculated critical points are denoted by the open circles. (The solid curves serve as a 
guide to the eye.) The dashed lines ((3e s — tu^/k^T — (3e s — e Ti )/k B T± 2S Byul ) correspond to 
the exact transition line of the closed packed trimer problem. 13 ' The dash-dotted lines correspond 
to the exact critical temperature Ae/k-sT c — ln[(ll + 5\/5)/2] 2.406 (Ae = 3e s — £h 3 , 3e s — £t 4 ) 
of the hard hexagon model. 18 ' . 



Clearly, the sublattice trimer concentration p takes a constant value p^ ~ 1/6 [ML] 
in the disordered phase. This constant behavior with p^ = 1/6 [ML] agrees with 
the exact solution. The chemical potential pLyN in Ref. Ill) is given by p,yN = 



(ct 4 — £Yt 3 )/kB,T. In Ref. |13|), it is shown that at py N = 2S sjm /kB, the system 
undergoes a discontinuous order-disorder phase transition. Hence, in our calculation, 
the transition temperature T m is given by -e T Jk B T m = (4/3) ln(3\/3/4) « 0.34883. 
The p-dependent free energy of Ref. H) has a cusp at p = pa, which suggests that 
the disordered phase persists as a metastable phase in a parameter region within 
the ordered phase near the transition line. In our calculation, the high-temperature 
fixed point of the PWFRG iteration, which corresponds to the disordered phase, 
actually continues to exist in the ordered phase to some extent. We thus see that 
our PWFRG calculation reproduces the exact solution. 



For general e s , the behavior of the system is governed by the ratio 

3e s — en 3 
r = -. 

3e s - e T4 



(3-1) 



3 ' 



The total trimer (adsorbate) concentration C tota i is given by C to tai = Ct 4 + Ch 
The monomer (free dangling bond) concentration is then given by 1 — 3C to tai- I n 
Fig. 6, we show the temperature dependence of Ct 4 , C*h 3 and C to tai f° r several 
values of r, obtained using PWFRG. This figure gives the relationship between the 
thermal adsorbate concentration and the adsorption energy. As a crosscheck, we also 
performed Monte Carlo calculations (not shown here). We find that the results of 
these calculation agree with those of the PWFRG calculation. 
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In our PWFRG calculation, we assume that the system is in the disordered 
phase, where the maximum eigenvalue eigenstate of the transfer matrix has transla- 
tional invariance. In the ordered phase.where the translational invariance is broken 
and the system has period-3 structure, Eif the PWFRG iteration ceases to converge. 
Hence, if the order-disorder transition is continuous, we can determine the transition 
line from the point at which the PWFRG iteration begins to lose convergence. By 
such a procedure, we obtain the phase diagram of the trimer-monomer system in 
Fig. 7. In the limit r—* -co or 1/r — > — oo the system can be described by the 
hard hexagon model, EJ) whose exact transition temperature is reproduced in our 
calculation. 

We should, however, note that to show whether the transition is continuous 
(second order) or discontinuous (first order) is a non-trivial problem. In the hard 
hexagon limit, the transition is known to be continuous,^ while, in the trimer-tiling 
limit (r — > 1), the transition seems to be first order. We do not discuss this point 
further, as it is beyond the scope of the present paper. However, we remark here 
that a Monte-Carlo calculation (not shown here) suggests that, at least for the cases 
in which r is not close to 1, the phase transition is second order, while as r — ► 1, 
the transition tends to exhibit a discontinuous (first-order- like) character. Detailed 
discussion of such points will be given elsewhere. 

§4. Application to Si(lll) lxl surface 

From the rocking curve of RHEED, Kohmoto and Ichimiya@) measured the 
adatom concentration and obtained Cr 4 = 0.20, Ch 3 = 0.05 and C to tai = 0.25 
at 900° C. Comparing these values with the results of our calculation, we determine 
the adsorption energies as 3e s — 6t 4 = 0.22 [eV/adatom] and 3e s — en 3 = 0.056 
[eV/adatom]. The difference en 3 — £t 4 = 0.17 eV is consistent with that estimated 
by Kohmoto and Ichimiya.LP 

We should note here that e s depends on the ambient vapor pressure through 
the surface chemical potential, as the Langmuir isotherm of monolayer adsorption, tf 
Hence ex 4 — 3e s and en 3 — 3e s will vary with the ambient vapor pressure, which may 
lead to a pressure- induced phase transition. 

§5. Summary 

In the present paper, we have discussed the generalized trimer tiling problem 
(trimer-monomer mixture problem) on a triangular lattice, where each trimer is 
interpreted as a triplet of dangling bonds satisfied by a trivalent adsorbed atom. 
This dangling-bond-based construction of the model, which may be regarded as 
an extension of the Sakamoto-Kanamori mappingllJ^ to the hard-hexagon model, 
makes it directly applicable to the Si(lll) lxl surface with both the H3 and T4 
sites considered as adsorption sites. We have determined the energy parameters of 
the system that consistently account for the results of experiments. Our scheme for 
the construction of the tiling model is also helpful for analysis of the model; we have 
introduced there exists a novel mapping to a square-lattice vertex model to which we 
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can apply the transfer-matrix method combined with the density-matrix algorithm. 

By considering adsorption of the n-valent adatom, we can extend the present 
analysis to more general problems, e.g. the "n-mer tiling" problem, the "n-mer and 
n'-mer mixture" problem, and so on. This extension is an interesting subject for fu- 
ture study. Clarification of the relation to other tiling problems^ 'EJ^ca^H'EfP'Ej)'^' 
is also an interesting topic to be explored. 
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